INTEGRAL REPRESENTATIONS FOR MULTIPLE 
,lj: ■ HERMITE AND MULTIPLE LAGUERRE POLYNOMIALS 

o 

S . PAVEL M. BLEHER AND ARNO B.J. KUIJLAARS 

C , Dedicated to Pierre van Moerheke on his sixtieth birthday 

f^ ' Abstract. We give integral representations for multiple Hermite and 

fT^ \ multiple Hermite polynomials of both type I and II. We also show how 

these are connected with double integral representations of certain ker- 
^ nels from random matrix theory. 

u 

j^ , 1. Multiple orthogonal polynomials 

Multiple orthogonal polynomials are an extension of orthogonal polyno- 
mials that play a role in the random matrix ensemble with an external source 

\0 . (1.1) ±e-Tr(y(M)-^Af)^^ 

\^ • defined on n x n Hermitian matrices, see |1J|S1E]. Here A is a fixed n x n 

^^ , Hermitian matrix and 1/ : M ^ M is a function with enough increase at ±00 

^T ■ so that the integral 

o ■ 

-(— > 

2 ! converges. Random matrices with external source were introduced and stud- 

S ■ ied by Brezin and Hikami [3 El El UHl E] , and P. Zinn- Justin [H CHI • 

^ . In what follows, we assume that A has m distinct eigenvalues ai, . . . , Um of 

L> I multiplicities rii, . . . , rim- We consider m fixed and use multi- index notation 

'j_| ■ n = (ni,. . . ,nm) and \n\ = ni -\ h rim- 

5t ! The average characteristic polynomial P^{x) = E [det(x/ — M)] of the 

ensemble (|l.ip is a monic polynomial of degree \n\ which satisfies for k = 
l,...,m, 

/oo 
Pfi{x)x^Wk{x)dx = 0, j = 0, . . . ,nfc - 1, 
-00 

where 

(1.3) u;fc(x) = e-(^(^)-"''-^'), 



Zr, = / e-T'-(^(^)-^^^)dM 
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see jl]. The relations p.2j) characterize the polynomial P^ uniquely. For 
^ = 0, we have that P^ is the usual orthogonal polynomial with respect to 
the weight e"'^^^''', which is a well-known fact from random matrix theory. 
For general m, the relations (|1.2() are multiple orthogonality relations with 
respect to the weights (|1.3() and the polynomial P^ is called a multiple 
orthogonal polynomial of type II. 

The multiple orthogonal polynomials of type I consist of a vector 

(1.4) {4\4\...,A^^\ deg4'^<n,-l, 

of polynomials such that the function 



(1.5) Qriix) = f24'\^)M^) 

fc=l 

satisfies 

(1.6) 1^ x^QH{x)dx = I J| ^: I °:- ^ '^1^1 - 2' 

The polynomials A\ are uniquely determined by the degree requirements 
1)1.4(1 and the type I orthogonal relations (|1.6|) . 

By the Weyl integration formula, the random matrix ensemble (|1.1() has 
the following joint eigenvalue distribution 

^ n ^'"^^"'^ ( I ^^^""^'^du] n(^^- - ^^)' d^^dX2 ■■■dXn 

^" i=l ^-^ ^ i<fc 

where dC/ is the normalized Haar measure on the unitary group U{n) and 
A = diag(Ai, . . . An). Using the confluent form of the Harish-Chandra / 
Itzykson-Zuber formula |13l I14j to evaluate the integral J e dU, we 

find that the joint eigenvalue distribution can be written as 

(1-7) -yr det (/i(Afc))i<j,fc<„ det (5j(Afc))i<j, fc<„ dXidX2 • • • dXn, 

where fj{x) = x^^^ and gi, . . . ,gn are the functions 

(1.8) 2;^e-(^(^)-'''=^\ k = l,...,m, j = 0,l,...,nk-l 

taken in some arbitrary order. Then ()1.7() is a biorthogonal ensemble in the 
sense of Borodin ^, which in particular implies that the eigenvalue point 
process is determinantal, that is, there is a kernel K{x, y) such that the k 
point correlation functions have determinantal form 

(1-9) det(K(A„A,)),<._^.<,. 

The multiple orthogonal polynomials of types II and I can be used to 
biorthogonalize the functions fj and Qj and to give an explicit formula for 
K. Choose a sequence of multi-indices no, ni, . . . , n^ = n such that \nj\ = j 
and fij < fij+i and define 

Pj (x) = Pn, (x) , Qj {y) = Qh^^^ . 
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Then pj+i G span{/i, . . . , /j}, and g^+i G spa,n{gi, . . . ,gj} for a suitable 
ordering of the functions (|1.8)) . In addition we have the biorthogonahty 



Pi{x)qj{x)dx = 5ij. 

o 

As in |S^ it then fohows that 

(1.10) K{x,y) = Y,pj{x)qj{y). 

j=0 

By the Christoffel-Darboux formula for multiple orthogonal polynomials ^ 
the kernel K satisfies 

(1.11) (x - y)K{x,y) = PH{x)Qn{y) - Y, -^Pn-s,{x)Qri+eM 

where 

(1.12) hf = / P^{x)x'''^Wk{x)dx 



and Cfc is the /cth standard basis vector in M™. 

In the following two sections we study two special cases related to multiple 
Hermite polynomials and multiple Laguerre polynomials. These cases corre- 
spond to the random matrix model ((171]) with V{M) = \M'^ and V{M) = M 
respectively (in the latter case we restrict to positive definite matrices). 
Brezin and Hikami ^8; and Baik, Ben Arous, and Peche ^ gave double in- 
tegral representations for the correlation kernels for these cases. We will 
derive integral representations for the multiple Hermite and multiple La- 
guerre polynomials, and use that to show that the kernels agree with the 
multiple orthogonal polynomial kernel 1)1.11(1 . 

2. Multiple Hermite polynomials 

The special case V{M) = ^M^ was considered in a series of papers of 
Brezin and Hikami, [Tj |H1 El llUllllj . This case corresponds to M = H + A 
where H is a random matrix from the GUE ensemble (1/Z„)e~2 dH 
and A is fixed as before. In |S] the following expression for the kernel was 
derived 
(2.1) 

^ ^^ {2mf J_,^ Jr l}^ \t-aj s-t 

where F is a closed contour encircling the points ai,...,am once in the 
positive direction, and the path from —ioo to ioo does not intersect T, see 
also Johansson |15j . 

When V{x) = 2^'^^ the multiple orthogonal polynomials are called mul- 
tiple Hermite polynomials, since they clearly generalize the usual Hermite 
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polynomials ^ [21 E|- We derive integral representations for the multi- 
ple Hermite polynomials of type I and type II, which resemble the integral 
representation (|2.H) of the kernel. 



2.1. Multiple Hermite polynomials of type II. The multiple Hermite 
polynomial P^ is the monic polynomial of degree \n\ that satisfies (|1.2|) with 
Wk{x) = e~2^ +akx^ k = 1,... ,m. 

Theorem 2.1. The multiple Hermite polynom,ials of type II has the integral 
representation 

(2.2) Pri{x) = -=: / e2(--)' J](s - auT" ds 

V2vrz J_ioo ^~^ 

Proof. Let us denote the left-hand side of (|2.2|) by P{x). After performing 
the change of variables s = t + x, we get 

(2.3) P{x) = ^ / 62*' TT (t + X - akT*^ dt, 

V27ri J-ioo ^J^ 

which shows that P is a polynomial of degree \n\ with leading coefficient 



1 ftOO 

62* dt = 1. 

ioo 



2m 



So P is a monic polynomial. 

Now we use (|2.3() to compute for /c = 1, . . . , m, and j = 0, 1, . . ., 

/oo 
P(x)x^e-2^'+«'=^dx = 
-oo 

(2.4) ^^ / / e3(*'-(^-"^)')x^JJ(t + x-a;)"'(ixdt. 

y ZTTI J —too J —oo 71 

Switching to polar coordinates x — a^ = r cos 6, t = ir sin 0, we find that the 

right-hand side of (|2.4() is equal to 

(2.5) 

r-27r 



eH f°° _i 2 



27r Jo 



/"CXD /"ZTT 

/ e-2'-V"^+i / (ak + r cos eye'"'^^T[ire'^ + ak-ai)'''de 



dr. 



The 0-integral vanishes for j = 0, ... ,71^ — 1, since the integrand can be 
written as a linear combination of e*^^ with integer p > n^ — j- Hence P is 
the multiple Hermite polynomial of type II and (|2.2() follows. D 
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Remark 2.2. Evaluating (|2.5j) for j = Uk we find that the ^-integral is 

2^(i)a^fcK-a;r and 



/oo 
-oo 



Uk^-^x^+akX^^ 



(2.6) = \/2^K)!eHj](a,-a,r- 

2.2. Multiple Hermite polynomials of type I. The multiple Hermite 
polynomials are polynomials A^ as in (|1.4|) such that the linear form 



(2.7) Qn{x) = Y,Af{x)e-^ 



fc=i 
satisfies (|1.6|) . 

Theorem 2.3. T/ie multiple Hermite polynom^ials of type I have the integral 
representation 

(2.8) ^?^(x)e-5^'+'^'=^ = -^ / e-5(*-^)' TT(t _ a^""' (it 

where T^ is a closed contour encircling a^ once in the positive direction, hut 
not enclosing any of the other points ai, I ^ k. 

In addition the function Qfi of 1)2. 7|) has the integral representation 

1 r "^ 

(2.9) Qn{x) = -j=—, i e-^(*"-)' n(* - «')""' dt 

where T is a closed contour encircling oi, . . . , Um once in the positive direc- 
tion. 

Proof. By the residue theorem, we have that the right-hand side of ()2.8() is 
equal to 



1 1 /" d 

(2.10) 



2TT{nk -1)! V'^* 






It is easy to see that (|2.1()j) has the form Ak{x)e~'i^ "'"'^'=^ where A^ is a 
polynomial of degree n^ — 1. Define the linear form 

(2.11) Qix) = V Ak{x)e-h-'+-''- = -=— I er^--)" [](i - a,)'^^ dt 

where T encloses all the points Oj, j = 1, . . . ,m, once in the positive direc- 
tion. 
Then 



1 /■ / 1 



oo 



x^Q(x)dx = — (^ ( -^ I " x^e-^^^'''^''dx]Y](t-ai)-'"dt. 
27^^ /r Vv2^7-oo J fJl 
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Since 

V 27r J-oo V 2vr J_oo 

is a monic polynomial of degree j in the variable t, we find 

/OO If ^^ 

-OO 27ri Jr ^^ 

where ttj is a monic polynomial of degree j. Deforming the contour T to 
infinity, and using the fact that the integrand is s-?"'"! + 0{s^~'"''^) as 
s — > OO, we find that 



/OO 
x^Q{x)dx = for j = 0, . . . , \n\ - 2 
-OO 



and 



/OO 
x^Q{x)dx = 1 for j = \n\ — 1. 
-OO 

This shows that Q = Qh and ^^ = A^}^ so that (|TH1) and (gill) follow. D 

2.3. The multiple Hermite kernel. Let us now show that the Brezin 
Hikami kernel (|2.1I) agrees with the multiple Hermite kernel (|l.llj) . To that 
end we compute ^ + ^ for the kernel (|2.1() in two ways. 

First we have 
dK dK 
dx dy 

1 f ^ I .1. i (s-x)2 - 1 ft-„)2 -TT /s-afc\ '^ -s + x + t-y 



/too /* / \ 

ds i dt e^(^--)^-i(*-^)^ n '-^ 
-ioo /r M; \t-akj 



= {x- y)K{x, y)-j:^ ds <f dt e^i^-^f-^i^^v)' J] ( F^ 

{2my J_i^ /r ^^^^ \t-ak 

The last double integral factors into a product of two single integrals, which 
by ((T2) and ((221) leads to 

(2.13) 1^ + 1^ = (a; - y)K(x, y) - Pri{x)QH{y). 

For the second way we evaluate ^ by noting that dxe2^^~''^' = —dse^'-^~^' 
and integrating by parts the s-integral 

dK 1 r°° _,_ / .. d _u.-.^-lU-.^ ^ fs-akV 1 



r ds i dt ^e'^i^-^^'-W-y? f\ ('_^^\ 
J^ioo Jt OS ;M V * - «fc / 



dx (27ri)2 y_.^ 7r ds ^Jl\t-akJ s-t 

-I />«oo /• 



(27ri)2 



ioo Jr 






m / N "fe 1 I "^ -I 

^ ' — _?^fc ]_ 

.t — ak J s — t \ ^-^ s — ttk s — t 

k=l Kk=l 
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Similarly, we use dye'^y^'y) = —Q^Q-2y^-y) ^ and apply integration by parts 
to the t-integral, to obtain 

dy {2mY J_i^ /r 



.t — Uk J s — t \ ^-^ t — ak s — t 



We add dTTIll and (tTT^ 



9^ + ^ = ,A./"""d.</dtei(-)M(-.) 



dx dy (27ri)2 7„j^ Jp 

(2.16) = -Vnfc— ^ / ds (^ dt e^(^-^)'-^(*-^)' 

f-^ (27rz)^ 7_i<^ /r 



fc=i 

^ni. — 1 

T—r / ,5 — uj \ - \s — uk 

X 



n 



For every A;, the double integral in (|2.16|) factors into a product of two single 
integrals, which by (|2.2I) and (|2.9|) are given in terms of multiple Hermite 
polynomials. It leads to 

dK dK *" 

^ fc=i 

From (fTT^ and (fTT7|) we get 

m 

{X - y)K{x,y) = Pn{x)Qn{y) -^nkPn-ek{x)Qn+ek{y)^ 

k=l 



which agrees with 1)1.11(1 since 

(2.18) ^^ = I(S7 

because of (|2.6p . 

3. Multiple Laguerre polynomials 
Complex Gaussian sample covariance matrices have a distribution 

(3.1) 1 g-Tr(E-iAf) (^g^ ^y ^M 

defined on n x ra positive definite Hermitian matrices M. The matrix M 
arises as M = XX^ where X is an n x (n +p), matrix whose independent 
columns are Gaussian distributed with covariance matrix S. Here p is a 
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non- negative integer. The distribution (|H.1|) is also called a Wishart ensem- 
ble. We assume that Tj~^ has eigenvalues f3i,...,/3rn > with respective 
multiplicities rii, . . . , rim- 

Writing S"^ = I - A, we see that (jHH} takes the form (frT|) with V{M) = 
M, but restricted to positive definite Hermitian matrices. It follows that 
the ensemble (|3.1|) can be described with multiple orthogonal polynomials, 
which in this case are multiple Laguerre polynomials |^ [21 1161 I17j . (To 
be precise, they are called multiple Laguerre II in [21 El to distinguish 
them from another generalization of Laguerre polynomials, called multiple 
Laguerre I.) 

The eigenvalues of M follow a determinantal point process on (0, oo) with 
kernel K(x,y) given by (jl.llj) 

m I.W 

(3.2) (x - y)K{x,y) = PH{x)Qri{y) - J^ -^PH-e,{x)Qri+eM 



ih 



(k) 

n-Sk 



where now P^ is the type II multiple Laguerre polynomial and Qfi{x) = 

^ A^ {x)x^e~^''{x) is the linear form involving the type I multiple La- 

fc=i 

(k) 

guerre polynomials A-. . 

Baik, Ben Arous and Peche ^ gave a double integral representation for 
the correlation kernel 

(3.3) A-(.,.) . plp£*/^*e"-.. (i)"'"'n (fff )"' 7^ 

where S and T are disjoint closed contours both oriented counterclockwise 
such that S encloses and lies in {s S C | 3f?s < min^ P^} and T encloses 
the points /3i, . . . , (3m and lies in the right half-plane. 

In view of 1)3. 3() and our experience with multiple Hermite polynomials 
we expect integral representations for the multiple Laguerre polynomials as 
well. We will see that this is indeed the case, and we use this to study the 
connection between the kernels 1)3. 2() and (|3.3() . It will turn out that the 
two kernels are equal up to a multiplicative factor x^y~'P. However, this 
difference does not affect the correlation functions (|1.9() . 

3.1. Multiple Laguerre polynomials of type II. The multiple Laguerre 
polynomial of type II is a monic polynomial Pf^ of degree \n\ such that 



(3.4) / P^(x)x^+Pe-^'=^dx = 0, 
Jo 



k = l,...,m, j = 0, ... ,nfc - 1. 



Theorem 3.1. The multiple Laguerre polynomial of type II has the integral 
representation 

(3.5) Pn{x) = ^.'|!)+f ^'' / e-.H"1-^-i n(^ - (^kT'ds 



MULTIPLE HERMITE AND LAGUERRE POLYNOMIALS 9 

where T, is a closed contour around oriented counterclockwise, which does 
not enclose any of the f3k 's- 

Proof. Denote the right-hand side of H3.5() by P{x) and write C = ^^. J-J" /_g -.n^. 
Then 

oo j „ m 

(3.6) P{x) = Cx-^Y,^ f s^"l"l"P"^ n^^ " l^kT'ds. 

AH terms in H3.6|) with j > |n| + 1 vanish by Cauchy's theorem, as weh as all 
terms with j < p — 1, as we can see by deforming the contour S to infinity. 
It follows that P is a polynomial of degree |n|, whose leading coefficient is 

Hence P is a monic polynomial of degree |n|. 

We now verify the orthogonality conditions (|3.4j) . Take k = 1, . . . ,m and 
let j = 0, . . . , n/c — 1. Then 

P{x)x^+Pe-^>'''dx = C (t ( x^e^'-^^^'^dx j s'^^^-P'^ JJ(s - f3i)'''ds 

where we have assumed that S is so that 3f?s < /3fc for every s G S. The 
inner integral is j!(/3fc — s)~"'~^i so that 

/ P{x)x^+'Pe-'^^''dx = Cj\ i s-l"l-P-i(s - hT''-^-^ Wis - PiT^ds 

The integrand in the last integral has no pole at /3fc, since j < n,fc — 1. At 
infinity the integrand behaves like s~^~p~'^. So by deforming the contour S 
to infinity, we conclude that the integral is zero. The conclusion is that P 
is the multiple Laguerre polynomial of type II and (|3.5|) follows. D 

3.2. Multiple Laguerre polynomials of type I. The multiple Laguerre 
polynomials of type I A\ , k = 1, . . . ,m, have degrees 

(3.7) deg4''^<nfc-l, k = l,...,m 
and are such that 





Qn{x) = 


m 

--Y^AfxPe-^-^ 

h 1 


satisfies 






(3.8) 


/>oo 

/ x^Qfi{x)dx = 
Jo 


r 0, j = 0, . . 

1 1, j = \n\ - 



1. 



There is also an integral representation for the multiple Laguerre polynomi- 
als of type I. 
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Theorem 3.2. The multiple Laguerre polynomials of type I have the integral 
representation 

(3.9) 4'^(-)-^«-^^^ = -iM^^4TT^ e-^H\^\^^-'ll{t-PO-^dt 

where r^ is a closed contour around (3^, which does not enclose nor any 
of the other points Pi, I ^ k. 

In addition the function Qfi has the integral representation 

(3.10) Qn{x) = -J^^4f^^^)^cfe-^H\^\+P-^f[{t-Piy-^dt 

where T is a closed contour around Pi, ■ ■ ■ ,Pm, but which does not enclose 
0. 

Proof. Only the pole at Pk contributes to the integral in (|3.9j) . By the residue 
theorem, we have that the right-hand side of (|3.9p 



„ . d 
const x^ , , 
dt 



rtfe-l 






which is easily seen to of the form Ak{x)x^e~^''^ where A^ is a polynomial 
of degree n^ — 1. 
Let 



fc=l 

which is then equal to the right-hand side of (|3.1flp . We may assume that T 



Q{x) = ^Ak{x)x'Pe- 

k=l 

le right-hand sid 
is entirely in the right half-plane. Then 

Jo ^ ' 27n{\n\+p-iy.Jr Jo l\^ 

The x-integral is (j + p)\t^^^'^^^, so that 



V "■■ 



1=1 



Assuming j < \n\ — 2, we can deform T to infinity without picking up a 
residue contribution at t = 0. The integrand behaves like t~^~'^ at infinity, 
and so ()3.1H) vanishes for j < \n\ — 2. For j = \n\ — 1, we pick up a residue 
contribution at t = 0, and the result is that 

(3.12) / xl"l-^Q(x)dx = 1. 

Jo 

Thus Q satisfies the type I multiple orthogonality conditions (|3.8|) and the 
theorem follows. D 
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3.3. The multiple Laguerre kernel. We finally compare the representa- 
tions (|3.2j) and (|3.3|) of the kernel. 

We start from the double integral (|3.3|) and evaluate xK{x, y) using an 
integration by parts on the s-integral. The result is 

1 /./..„..-.* ^A'"'^" A ^-A^"' 






J -Pi 



1 I Inl + n v-^ nu 1 1 

Similarly, after integration by parts on the i-integral, 

1 \ \n\ + p ^^ nu 1 I 
3.14 X <^ ^ - > V + } 

s-t\ t ^i-A s-tj 

Hence 

{x -y)K{x,y) 

I \n\ +p , sr^ rik 



St 



'^ti'i^-Mit-Pk) 



n\+p-l ™ /„ _ /3, \ "' 



t'f/^--"^n{^) 



nk I ,I,.^s-,t f t\\'^^' rj f s - PiY^ {s-P,r- 



Eia,/.-//--"© n(^)' 



Now we have m + l double integrals and they all factor into products of two 

single integrals of the forms 1)3. 5|1 and 1)3. 1U() . The result is that 

(3.15) 

{X - y)K{x,y) = X^y'^ Pn{x)QH{y) - ^^k ^ Pn-ekix)Qn+ekiy) ■ 

V fc=l ^k J 



It can be shown that 



\n\+p hf 
rik- - 



Pi hi'K 

SO that (|3.15j) agrees with (|3.2j) up to the factor x^y~'P. However, this factor 
is not essential since it does not change the correlation functions (|1.9j) . Hence 
(|3.2j) and (|3.3|) are essentially the same. 
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